Dynamics of Few Co-rotating Vortices in Bose-Einstein Condensates 
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We study the dynamics of small vortex clusters with few (2-4) co-rotating vortices in Bose-Einstein 
condensates by means of experiments, numerical computations, and theoretical analysis. All of 
these approaches corroborate the counter-intuitive presence of a dynamical instability of symmetric 
vortex configurations. The instability arises as a pitchfork bifurcation at sufficiently large values of 
the angular momentum that induces the emergence and stabilization of asymmetric rotating vortex 
configurations. The latter are quantified in the theoretical model and observed in the experiments. 
The dynamics is explored both for the integrable two- vortex system, where a reduction of the phase 
space of the system provides valuable insight, as well as for the non- integrable three- (or more) 
vortex case, which additionally admits the possibility of chaotic trajectories. 



Introduction. The realm of atomic Bose-Einstein con- 
densates (BECs) [H has offered a pristine setting for stud- 
ies on the dynamics of few- vortex clusters f2| . Most inves- 
tigations, however, have focused on either a single vortex 
or large scale vortex lattices Recently, theoretical 

investigations on the study of clusters of 2-4 vortices 
15[, have appeared, chiefly motivated by the experimen- 
tal realizations of such states 16 - 1^ . This focus has been 
heretofore centered on the fundamental building block of 
the vortex dipole, i.e., a pair of counter-circulating vor- 
tices. 

Our aim in the present work, in contrast, is to ex- 
plore the dynamics of small vortex clusters of 2-4 vor- 
tices that belong to the co-rotating (same charge) va- 
riety. The original work of Ref. 0| and subsequent ef- 
forts [20i] have already paved the way for an understand- 
ing of symmetric few- vortex configurations rotating as a 
rigid body, and their three-dimensional generalizations, 
i.e., U- and S-shaped vortices, as well as vortex rings 21 1. 
In this context, our work presents a rather unexpected 
twist: we have found that, under suitable conditions, the 
usual symmetric, co-rotating vortex configurations (cen- 
tered line, triangle, and square) become dynamically un- 
stable. More specifically, these states become subject to 
symmetry breaking, pitchfork bifurcations that lead to 
the spontaneous emergence of stable asymmetric rotat- 
ing vortex clusters. 

We present our analysis of these features in the in- 
tegrable (at the reduced particle level) setting of a co- 
rotating vortex pair, and illustrate their generality by 
further considering a rigidly rotating vortex triplet and 
quadruplet. In the first case, we devise a theoretical for- 
mulation that not only explores the instability and man- 
ifests its growth rate, but also enables a visualization 
of a two-dimensional reduced phase space of the system 
in which the pitchfork bifurcation becomes transparent. 



In the latter cases, we suitably parametrize the system, 
exploring the different regimes of symmetric and asym- 
metric periodic orbits. Our theoretical analysis treats 
vortices as classical particles, with dynamics governed by 
ordinary differential equations (ODEs). This reduction of 
the original vortex cluster system allows for the analytical 
characterization, numerical observation and experimen- 
tal confirmation of the symmetry breaking phenomena. 



Theoretical Analysis. As illustrated in Refs. [15 
IToj , and justified by means of a variational approx- 
imation |22l |. vortex dynamics governed by the two- 
dimensional mean-field Gross-Pitaevskii equation. 
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can be reduced to a system of ODEs for the vortex posi- 
tions. In the original partial differential equation (PDE) 
model (jlj, the time is measured in units of while 
the positions are measured in units of harmonic oscilla- 
tor length along the z-direction and fl = ujx/tJz = tUy/uz, 
with ujj being the harmonic trap frequency along the j- 
direction (see, e.g., Ref. [4j). This ODE reduction is the 
starting point for our analysis in the co-rotating case. 

The dynamics of vortex m at position {xm,ym) arises 
from two contributions: (i) a position-dependent vortex 
precession about the trap center with frequency Sm^pr, 
and (ii) a vortex-vortex interaction with vortex n that 
induces a velocity perpendicular to their line of sight of 
magnitude S'„ Wvort/Pmrn where p„i„ is the distance be- 
tween vortices m and n, Sm and Sn are their respec- 
tive charges, and Wvort is a dimensionless constant; see 
Ref. la, ll9| . The equations governing the dynamics of N 
interacting vortices embedded in a condensate are there- 
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The precession about the trap center can be approxi- 
mated by Wpr = Wpi./(1 — ^^/^tf)' where the frequency 
at the trap center is cUp^. — In (A^) /i?|.p, /i is the chem- 
ical potential, Rtf — \/2jI/n is the Thomas- Fermi (TF) 
radius, and A is a numerical constant I, [H, [H. To 
describe better the actual vortex dynamics in the trap, 
the constant Wvort in Eqs. ^ may be adjusted to ac- 
count for the screening of vortex interactions due to the 
background density modulation [23 1. 

We now focus on the case of two identical vortices of 
unit charge Si = S2 = We proceed to adimensionalize 
Eqs. ([2]) by scaling {x,y) by Rtf and time by l/i^pr, and 
use polar coordinates {xn,yn) = (''n cos(0„), r„ sin(0„)). 
We then seek symmetric stationary states ri = r2 = 
and Oi — 02 ~ TT, and find the following frequency of the 
co-rotating vortices: 
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where c = |(wvort/Wpj.) yields a measure of the rela- 
tive strength of vortex interaction and spatial inhomo- 
geneity. The comparison of the orbital frequency be- 
tween the ODE and the PDF models is given in Fig. [T^. 
Given the co-rotating nature of this state, consideration 
of 5mn ~ 0m ~ renders this state a stationary one; lin- 
earizing around it using r,„ = r, 4- Km and dmn — + Sm 
yields the following equations of motion for the pertur- 
bations about the symmetric equilibrium: 
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7 — ^rnr ■ It is then straightforward to ob- 

2r^ (l-r^) ^ 

serve that this squared epitrochoidal (motion of a point 
in a circle that is rotating about another circle) rela- 
tive precession frequency of the two vortex positions and 
phases changes sign at r^j. — sfcji^sfc^ 2). This signals 
our first fundamental result, namely the destabilization of 
the symmetric 2-co-rotating vortex state for sufficiently 
large symmetric distances of the vortices from the trap 
center. A comparison of the ODE and PDF models for 
the orbital and epitrochoidal precession frequencies for 
these two cases is given in Fig. [T^-b. 

The dynamical instability of symmetric states suggests 
the potential existence of additional, asymmetric, ones. 
Seeking states with dmn — tt and r* ^ yields 

-rlr;{rl+r;r + c{l-rl') (l - ) = 0, 



which will be the condition defining our radially asym- 
metric co-rotating solutions. The mirror symmetry of 
the 2-vortex system predisposes towards the pitchfork, 
symmetry breaking nature of the relevant bifurcation, a 
feature verified by the diagram of Fig. [TJ;. This diagram 
is given for the angle (f) = tan^^ r2/ri as a function of the 



angular momentum Lq = rf 



'2: 



which is a conserved 



quantity for our system. It is interesting to note that if 
the single dimensionless parameter of the system is small 
(c < 3), then the critical value Lcr for Lq — at which 
the bifurcation from symmetric to asymmetric periodic 
orbits occurs — is supercritical, while if c is sufficiently 
large (c > 3), it becomes subcritical [11] (not shown). 

To elucidate the pitchfork nature of the bifurcation, we 
develop a phase plane representation for all 2-vortex con- 
figurations. The integrability of the reduced 2-particle 
description can be understood on the basis of the fact 
that this 4-dimensional system has two integrals of mo- 
tion, namely the angular momentum Lq, defined above, 
and the Hamiltonian H, which can be written in polar 
coordinates as 



1 



H = -ln[{l-rl) (1-ri)] - -In [r? 



where D = 2r\r2 cos((5) and b — B2~Q\. Using Lq and the 
angle to express r\ and r2, one can rewrite the Hamil- 
tonian as a function of ((/), ^) having thus effectively re- 
duced the 4-dimensional system into a 2-dimensional one. 
Thus, for different values of Lqj we can represent the or- 
bits in the effective phase plane of (</>, S) in which the dif- 
ferent orbits correspond to iso-energetic if((/), (5)=const. 
contours. This is done in Fig. [5] for values that are both 
below and above than the critical value of Lq at fixed c. 
It can then be inferred that the symmetric fixed point 
with ((/), (5) = (7r/4,7r) is stable in the former case, while 
it destabilizes in the latter case through the emergence 
of two additional asymmetric (^ 7^ '^1^ states along the 
horizontal line J = vr of anti-diametric vortex states. 

Remarkably, although the properties of the system dra- 
matically change as we go from two vortices to three and 
four, the symmetry breaking bifurcation associated with 
the symmetric solutions persists. In particular, when 
TV > 2, the persistence of the two conservation laws 
discussed above is not sufficient to ensure integrability 
of the system, and its absence is manifested in a dra- 
matic form in the resulting 6 (iV = 3) and 8 (iV = 4) 
dimensional systems through the presence of chaotic or- 
bits. Nevertheless, one can still theoretically analyze the 
highly symmetric co-rotating states of the system. 

For = 3, this state is an equilateral triangle such 
that r\ — r2 — Tj, = r* and (J^.i+i = 27r/3, with an or- 
bital frequency predicted as aJorb,3 = + jztt- In the co- 
rotating frame, the linear stability analysis around this 
rigidly rotating triangle can be performed giving rise to 
an epitrochoidal frequency w^p g = — {i-r^y^ • -'^^ ^'^^^ 
case too, a critical radius exists r^j. 3 
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FIG. 1: (Color online) (a) Orbital and (b) epitrochoidal frequency as a function of the radial position from the trap's center 
for two vortices. Both frequency and radial position are in rescaled units. The solid line represents results from the ODE and 
the dotted line from the PDE. The vanishing of the latter signals the onset of instability. Here, — 0.05 and /i = 1. Panel 
(c) depicts the quantity (/f>/7r, which equals 1/4 when ri = r2, as a function of the square root of the angular momentum for 
c — 0.1. Panels (d) and (e) depict the corresponding phenomena for A'' = 3 and A*' = 4 vortices for c — 0.1. Panels (c), (d) 
and (e) include a few configurations along the main bifurcation branches ([blue] solid and [red] dashed lines corresponding, 
respectively, to stable and unstable configurations) depicting the relative position of the vortices (red triangles) with respect to 
the center of the condensate (green crosses) . 




FIG. 2: (Color online) (a) Contours from the reduced Hamil- 
tonian for vortices close to the center of the trap, i.e. Lo = 0.16 
below critical (Lo < Lcr ~ 0.273), and (b) their corresponding 
orbits. Panels (c) and (d) depict, respectively, the same quan- 
tities but for vortices further out from the center for Lo — 0.36 
above critical (Lo > Lcr = 0.273). The circles and triangles 
correspond, respectively, to the initial and final positions of 
the two vortices. Each color on the contour plot matches the 
corresponding orbit in the position diagram. Here, c = 0.1. 



such that the symmetric state is destabilized and asym- 
metric orbits arise and are stable past this critical point 
as can be seen in Fig. [TJi. The dynamical picture is con- 
siderably more complicated but the conservation of the 
angular momentum ensures that the dynamical evolution 
resides on the surface of a Bloch sphere. We thus define 
two angular variables tan0 = r2/ri and cos 6* — r^/^/Lo 
and depict the associated pitchfork bifurcation in Fig.[TJi 
for the subspace of solutions constrained to ri = r2 and 
^12 = ^23- This bifurcation diagram describes a vor- 
tex configuration containing a stable symmetric rotat- 
ing triangle before the bifurcation and stable asymmet- 
ric rotating triangles after the bifurcation. In addition to 
the equilibrium and near-equilibrium orbits, we observe 



chaotic orbits arising both in a more localized form, ex- 
ploring the vicinity of equilibrium orbits, and in a more 
extended one spanning all space (not shown). 

While the general phenomena for iV = 4 are already 
rather complex, some basic features can still be inferred 
and the symmetry breaking nature of the proposed insta- 
bility persists — cf. Fig. [TJs. Here, (j) = tan~^r3/ri, and 
we have constrained the vortices to be in a cross with 
right angles and ri = and r2 = r^. A general ex- 
pression for the orbital frequency of the rigidly rotating 
state is uJorh,N = ''^"j"'*'' -I- which is valid for any 

N. In the case of the square configuration with r.i — r^, 
and = 7r/2, there emerge two epitrochoidal vibra- 

tional motions with frequencies \/—Xi and \/—X2, where 
Ai = - J^, and A2 = - ^. These, in 



turn, correspond to two critical points: one identical to 
the one given above for the iV = 3 case, and one that 
is always higher, given by r^j. 4 — VSc/ {V^c + 2); hence, 
the same phenomenology persists. 

Experimental Observations. We now briefiy discuss 
experimental manifestations of the symmetry breaking 
events discussed above and of the emergence of asym- 
metric configurations. 

The details of the experimental setup may be found 
elsewhere 17 , 1^ . We begin with a magnetically-trapped 
BEC of TV ~ 5-8 X 10^ atoms in the |F = 1, mp = -1) 
hyperfine level of ^^Rb. The radial and axial trap fre- 
quencies are {u!r,i^z)/2TT = (35.8, 101.2) Hz. Vortices are 
introduced through a process of elliptical magnetic trap 
distortion and rotation [23] during evaporation [2^. In 
terms of the trap frequencies along the major and minor 
axes of the distorted potential, uj^ and ujy respectively, an 
ellipticity e = {lo^ — w^)/(w^ + w^) = 0.20 and a rotation 
frequency of 8.5 Hz usually produces a co-rotating pair. 
Higher rotation frequencies are used to generate larger 
numbers of co-circulating vortices. 

A partial-transfer (5%) imaging method 
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ployed to create a sequence of atomic density profiles. 
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as shown in Fig. [3^-d. The effect of the extractions is 
primarily to diminish the number of atoms in the con- 
densate 3 12 1- Curiously, atomic losses have little ef- 
fect on the calculated coupling parameter c, which scales 
only as log TV; thus c falls between 0.11 and 0.10 over the 
range N = 0.3-0.8 x 10^ atoms. For convenience, we take 
c = 0.1 in the following analysis. 

We examine 52 experimental time series, each consist- 
ing of 8 snapshots spanning a total time of 240 to 480 ms. 
For each snapshot the center of the vortices and the ra- 
dius of the BEC cloud are extracted using a least square 
fitting algorithm [see panels a)-d) in Fig. [3]. The vortex 
positions are then normalized to the BEC radius (i.e., 
TF units) and the angular momentum Lq and Hamilto- 
nian H were computed for each frame [see panels e)-h) 
in Fig. [3]. For each series, the averaged angular momen- 
tum Lq and Hamiltonian H are computed (see horizontal 
dashed lines in the middle panels in Fig.|3]). Using Lq we 
compare the experimental points representing each orbit 
in the {(j), S) plane to the isocontour of H corresponding 
to H, as shown in the right column of panels in Fig. [31 
and find good agreement between the two. 

The panels a)-d) in Fig. |3] depict a cohort of typical 
time series, together with their respective fits, that ex- 
emplify the different qualitative cases that we observed 
in the experiments. In particular, we find that the dy- 
namics of the vortices depends on whether the average 
angular momentum is below or above the critical thresh- 
old Lcr = 2r^^. This separates cases where asymmetric 
orbits are, respectively, non-existent and possible. The 
different qualitative cases that we observe, which are is 
displayed in Fig. [31 may be grouped as follows: 

• For Lq < Lcr and relatively small H, the experiment 
displays symmetric orbits. See rows a) and e) in Fig. [31 

• For Lq > Lcr and moderate H, the experiment dis- 
plays (i) symmetric orbits where both vortices (in the 
co-rotating frame) are approximately on the same side of 
the cloud chasing each other on the same path [see rows 
b) and f ) in Fig. [3] or (ii) asymmetric orbits [see rows c) 
and g) in Fig. [3]. The choice between these two orbits 
is determined by the initial conditions. Initial conditions 
inside the area delimited by the separatrix (red double- 
loop curve in the right panels of Figs. [31 f and Figs. [31 g) 
emanating from the saddle point {(f), 6) — (7r/4,7r) give 
rise to asymmetric orbits. 

• For Lq > Lcr and large H, the experiment displays or- 
bits in which one vortex remains close to the center while 
the other orbits around it close to the periphery of the 
cloud. See rows d) and h) in Fig. [31 

As is clear from these examples and the remaining 48 
data sets that we studied (see supplemental material), 
asymmetric orbits are only found when Lq > Lcr and 
when the vortex orbits fall inside the asymmetric min- 
ima regions of the Hamiltonian picture in the {(f>, 6) plane. 
Asymmetric solutions absent in all of the cases for which 
Lq < Lcr- These results are in good agreement with 
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FIG. 3: (Color online) a)-d) Typical experimental series for 
the dynamics of two co-rotating vortices (time indicated in 
ms). Tiie large (green) circles and the (red) crosses represent, 
respectively, the fitted TF radius and center of the cloud while 
the small (yellow) dots depict the fitted vortex centers. The 
(red) dashed circles represent the critical radius, rcr, above 
which symmetric orbits become unstable. e)-h) Manifesta- 
tion of the pitchfork bifurcation for the experimental series 
depicted in panels a)-d) which correspond to c = 0.1. Left 
column: experimental vortex positions and their correspond- 
ing orbit from the reduced ODE model (solid line), in TF 
units in the co-rotating frame. Middle column: correspond- 
ing Lo (blue circles) and —H (green squares) and their av- 
erages (horizontal dashed lines) as well as the critical value 
for Lo (solid horizontal line). Right column: corresponding 
orbits in the {(f), 5) plane along with isocontours for constant 
H (highlighted in dark gray is the isocontour corresponding 
to the average H and in red is the separatrix delimiting the 
area containing asymmetric orbits). 



the theoretical prediction of the pitchfork bifurcation de- 
picted in Fig. [TJ;. 

To extend our considerations, we briefly present a com- 
parison between experiment and theory for = 3 and 
iV = 4 vortices. The main phenomenology is depicted 
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FIG. 4: (Color online) Experimental series for trios [a) and b)] 
and quartets [c) and d)] below [a) and c)] and above [b) and 
d)] the critical threshold. e)-h) Corresponding time series for 
I/O and —H. Same notation and units as in Fig. [3] 



using two examples for each case in Fig.U] Panels a) and 
b) correspond to the = 3 vortex case below and above 
the pitchfork bifurcation [see panel e) and f)]. Panels c) 
and d) depict the equivalent scenario for = 4 vortices. 
As the figure illustrates, and is observed in all of the 
cases that we studied (17 data sets for = 3 and 5 data 
sets for A'^ = 4; not shown), the main phenomenology for 
A^ — 3 and A^ = 4 persists in that all configurations with 
Lq > Lcr are not symmetric and symmetric configura- 
tions, or epitrochoidal oscillations about them, are only 
present when Lq < Let- 

Conclusions. We have revisited the theme of co- 
rotating few vortex clusters in atomic Bose-Einstein con- 
densates. By a combination of theoretical analysis, 
numerical computation and experimental observation, 
we have illustrated a strong manifestation of symmetry 
breaking through a pitchfork bifurcation, which led to 
the destabilization of symmetric solidly rotating config- 
urations and gave rise to the the emergence of stable 
co-rotating but asymmetric vortex configurations. We 
showed that this analysis is fruitful not only for the inte- 
grable (at the reduced particle level) two-vortex setting, 
where a suitable parametrization of the phase space was 
provided, but also for the non-integrable cases of A^ = 3 
and A^ = 4 vortices where chaotic orbits exist. 

Naturally, it would be interesting to provide a more 
global characterization of the dynamics of the three-body 
problem, which is perhaps the most analytically tractable 
and interesting case due to its potential for chaos. An- 
other expansion of the present considerations involves 
their generalization to higher dimensions. In this case, 
it would be interesting to see, upon gradual decrease of 
the trapping frequency in the third dimension, whether 
the symmetry-breaking phenomena persist for line vor- 
tices and vortex rings. These aspects are presently under 
study and results will be reported elsewhere. 
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-0.0125 
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0.282 


-0.111 
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-0.0216 
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-0.105 
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-0.0177 
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0.300 


-0.075 


1) 
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-0.0162 
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-0.137 


m) 
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-0.0294 
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-0.135 
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-0.0230 
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0.333 


-0.157 
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oo) 
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-0.166 
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-0.0273 
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-0.185 


q) 


0.179 


-0.0378 
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-0.196 
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0.184 


-0.0404 


rr) 
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-0.241 
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0.187 


-0.0429 
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0.480 


-0.286 
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0.193 


-0.0392 


tt) 


0.519 


-0.335 


u) 


0.194 


-0.0517 
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0.558 


-0.372 


v) 


0.194 


-0.0509 


vv) 


0.619 


-0.382 


w) 


0.205 


-0.0617 


ww) 


0.655 


-0.471 


x) 


0.208 


-0.0634 


xx) 


0.664 


-0.489 


y) 


0.213 


-0.0488 


yy) 


0.690 


-0.606 




0.213 


-0.0697 


zz) 


0.751 


-0.478 



TABLE I: Averaged angular momentum La and averaged 
Hamiltonian H for all the experimental series containing two 
co-rotating vortices. The critical angular momentum for the 
experimental setup is Lcr = 0.2731. 
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SUPPLEMENTAL MATERIAL: 

Here we present the results that we obtained for the 
52 experimental series containing two co-rotating (same 
charge) vortices. The experiments are ordered by the av- 
erage angular momentum Lq- which is computed from 
the fitted centers and cloud diameter. The average is 
then computed by averaging over the eight experimental 
snapshots for each series. TableUlists the average angular 
momentum and average Hamiltonian for all the experi- 
mental series containing two co-rotating vortices. The 
corresponding snapshots for all the series are shown in 
Fig.[5]using the same ordering as given in TableU Finally, 
Fig. [5] depicts the analysis for each individual series by 
presenting the orbits in the co-rotating frame (respective 
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FIG. 5: Experimental time series corresponding to the 52 experiments containing two co-rotating vortices. The snapshots 
were taken at regular times for a total time spanning 240 ms or 480 ms (for individual times for each snapshot refer to the 
corresponding middle column in Fig. |6]). The experiments are ordered according to average angular momentum as listed in 
Table in For ease of presentation, all the images have been centered and rescaled so that they occupy approximately the same 
visual area. 



left columns), the times series for the angular momen- 
tum and Hamiltonian (respective middle columns), and 
the orbits in the reduced (0, 5) plane (respective right 
columns). 



As observed in Fig. (SJ the experimental orbits in the 
co-rotating frame (respective left columns) match well 
the theoretical orbits from our model (see gray orbit) for 
all of the experimental runs. Significantly also, all ex- 
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FIG. 6: (Color online) Analysis for experimental series depicted in Fig.O For each series the respective columns represent: (i) 
The left column depicts the vortex orbits (green squares and blue circles), and their corresponding orbit from the reduced ODE 
model (gray solid line), in TF units in the co-rotating frame, (ii) The middle column depicts the time series for the angular 
momentum Lq (blue circles) and negative Hamiltonian —H (green squares) and their respective averages (horizontal dashed 
lines) as well as the critical value for Lo (solid horizontal line), (iii) The right column depicts the orbits in the {4>,S) plane 
along with the isocontours for constant H (highlighted is the isocontour corresponding to the average H). 



perimental orbits in the reduced (0, S) plane (respective 
right columns) match well the theoretical orbits corre- 
sponding to the isocontour of constant H equal to the 
average H from the experiment — with the exception of 
the experiment^rr which we discuss below. In most of 



the experimental series both the angular momentum and 
the Hamiltonian are approximately conserved, buttress- 
ing the validity of our dynamical reduction. 

To further bolster the case for the validity of our re- 
duced dynamics we now concentrate on the main result 
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of our manuscript. The reduced model predicts that, 
as the angular momentum is increased, the co-rotating 
symmetric configuration loses stability in a pitchfork, 
symmetry-breaking bifurcation. At the bifurcation point 
a pair of stable, asymmetric co-rotating configurations 
emerge. Here, the bifurcation point is predicted to oc- 
cur at Lcr = 2r^j. where the square of the critical radius 
is = \fcl{\fc -I- 2) in units of the Thomas- Fermi ra- 
dius. In our experiment, the value of c, the ratio be- 
tween the rotation induced by vortex-vortex interactions 
and the rotation induced by the precession of vortices 
due to the inhomogeneity of the cloud's background pro- 
duced by the external magnetic trap, is c w 0.1. This 
gives a value for the critical angular momentum in our 
experiments of = 0.2731. Therefore, corroboration of 
our analytical prediction about the pitchfork bifurcation 
would imply that (i) for all the experimental series with 
an average angular momentum below critical (Lq < Ld ) 
there should be no co-rotating asymmetric orbits, and (ii) 
if an asymmetric co-rotating orbit is observed its corre- 
sponding angular momentum should be above the critical 
threshold. The critical threshold for the angular momen- 
tum lies between experiments#hh and ii. 

It should be noted at this stage that an angular mo- 
mentum above critical does not imply an asymmetric 
co-rotating orbit. For example orbits with high angular 
momentum (see experiments#ss-zz) might correspond to 
(i) orbits where vortices seem decoupled and one vortex 
remains close to the trap center while the other orbits 
around close to the periphery of the cloud (see experi- 
mcnts#ss, tt, uu, ww, xx, yy), or (ii) the two vortices 
belong to the same path that has them on the same side 
of the cloud (see experiments#vv and zz). The type of 
resulting orbit is heavily dependent on the initial condi- 
tions. In fact, in the {4>,6) plane, the area of the set of 
initial conditions that leads to co-rotating asymmetric or- 
bits (or slight perturbations thereof) is rather small and 
therefore asymmetric co-rotating instances are less com- 
mon. Nonetheless, we were able to identify two orbits 
as being asymmetric (i.e., epitrochoidal perturbations of 
the asymmetric orbit) in experiments#oo and qq. In 
fact experiment^oo is featured in panels c) and g) in 
Fig. 3 of the manuscript. It is reassuring that the angular 
momenta corresponding to the experimental asymmetric 
cases is clearly above critical supporting our dynamical 
reduction and the pitchfork bifurcation picture that it 
unveiled. 

It is important to note that one of the 52 experimental 
series did not match our theoretical expectations. This 
case is experiment #rr. As it is clear from the right 
panel of Fig. |6lrr, the experimental orbit lies close to 
the point (0, <5) = (7r/4, tt) indicating an almost per- 
fect symmetric orbit as observed in the left panel of 
Fig. [6lrr. However, this case has an average angular mo- 
mentum of Lq = 0.438, well above the critical threshold 
Lcr = 0.2731. The reduced dynamics description predicts 



that this orbit should be unstable, lying as it does close to 
an unstable saddle. It is not clear how this symmetric or- 
bit can exist in this supposedly unstable region. It is pos- 
sible, albeit unlikely, that this is an unstable orbit that 
has not had enough time for the instability to develop. 
A second, more plausible, explanation can be given if 
one carefully examines the experimental snapshots cor- 
responding to this experiment depicted in Fig. [5lrr. It 
is evident that in this experimental series, the traces of 
the vortices are much more diffuse and blurry than in 
the other series. The vortex lines corresponding to the 
vortices in this series may not be straight or aligned with 
the trap axis but rather bent 2^ 27 1. This explains on 
one hand the blurriness of the vortex cores in the ex- 
perimental snapshots and the unusual (with respect to 
our analytical dynamical reduction) dynamics that they 
display. This would further serve to illustrate the conse- 
quence of the true three-dimensional nature of the con- 
densate that, in some instances, cannot be captured by 
our reduced two-dimensional description. 
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